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Abstract 

A theory of electron spin relaxation in semiconducting carbon nanotubes is developed based on 
the hyperfine interaction with disordered nuclei spins 1=1/2 of ^^C isotopes. It is shown that strong 
radial confinement of electrons enhances the electron-nuclear overlap and subsequently electron spin 
relaxation (via the hyperfine interaction) in the carbon nanotubes. The analysis also reveals an 
unusual temperature dependence of longitudinal (spin-flip) and transversal (dephasing) relaxation 
times: the relaxation becomes weaker with the increasing temperature as a consequence of the 
particularities in the electron density of states inherent in one-dimensional structures. Numerical 
estimations indicate relatively high efficiency of this relaxation mechanism compared to the similar 
processes in bulk diamond. However, the anticipated spin relaxation time of the order of 1 s in 
CNTs is still much longer than those found in conventional semiconductor structures. 

PACS numbers: 85.35.Kt, 85.75.-d 
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Due to their unique electrical properties, carbon nanotubes (CNTs) are considered to 
be the ultimate structure for continued "scaling" beyond the end of the semiconductor mi- 
croelectronics roadmap.— Moreover, the unique electrical properties of CNTs are enhanced 
by the equally unique structural properties. This combination assures the development of 
CNTs for important applications and has largely been the focus of attention to date (see 
Refs. y and |3| as well as the references therein). Recently, however, the researchers are 
beginning to explore other important advantages that the CNTs can offer. For example, 
the CNTs with naturally low or no impurity incorporation allow, in addition to the more 
conventional scaled transistor application, the injection and use of electrons with polarized 
spin^ as an added variable for computation. Thus, CNTs are an ideal medium for the 
development of the emerging field of spintronics.— Further, the anticipated long spin relax- 
ation times allow coherent manipulation of electron spin states at an elevated temperature, 
opening a significant opportunity for spin-based quantum information processing. Clearly, 
spin dependent properties of CNTs warrant a comprehensive investigation from the point of 
view of fundamental physics (see, for example, Refs. andO) and practical applications. 

The objective of the present paper is to theoretically investigate the electron spin relax- 
ation properties in the CNTs, a crucial piece of information for any spin related phenomena. 
Specifically, we consider the electron hyperfine interaction (HFI) with nuclear spins / = 1/2 
of ^^C isotopes (with the natural abundance of 1.10%). The HFI is thought to be one of the 
most important spin relaxation processes in the CNTs; strong radial confinement of electrons 
in the CNTs enhances electron-nuclear overlap and subsequently the hyperfine interaction 
compared to the bulk crystals. On the other hand, the mechanisms related to spin-orbital 
interaction are expected to be extremely weak in CNT.!'-^ In the following analysis, the main 
emphasis will be on the single-walled semiconducting nanotubes. 

The property of our interest is the longitudinal (Ti) and transversal (T2) spin relaxation 
time of an electron with the radius vector 7^ and spin ^ in a CNT. The governing Hamil- 
tonian caused by the Fermi contact HFI with nuclear spins / located at lattice sites 
can be expressed as 

N 

Hhf = Qoahf 5^7'5(r - it,) = S^, (1) 

where the HFI constant ah/ and the area of the graphene sheet Qo are normalized per carbon 
atom. As indicated, this Hamiltonian Hhf can also be expressed in terms of the fluctuating 
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field operator B that mediates spin relaxation. 

To proceed further, G must be expressed in terms of electronic Bloch states of the relevant 
energy bands. In an effective mass approximation, the eigenstates for the conduction bands 
in the vicinity of the K point take the forn>2*ii 

^ Tib ("f* 

where Aq denotes the length of the CNT, L (= ni 61 + 722 62) the chiral vector in 
terms of primitive translation vectors 61, h 2 and integers ni, ^2, and B^^n{k) = [^u,n — 
ik]/ ■>/ ae^^n^ + k"^ with se^n = {27r/L){n — z^/3); the quantum number n = 0, ±1, ±2, ... dis- 
tinguishes the energy bands, while u takes one of the three integers —1,0,1 that makes 
{rii —n2 — i') an integer multiple of 3. As shown in Fig. 1, ^ and r] represent the coordinates 
for the axes directed along L (i.e., the circumference) and the CNT (i.e., il?), respectively. 
The eigenstates \k') for the K' valley can be readily obtained from Eq. (j21) by substituting 
B*_u,w ^'^d k — > The wave vectors k and k' are determined from the K and K' 
points of the Brillouin zone, respectively. 

The corresponding dispersion relation for the \k) states reads 



£n,k = 7V'^^a™^"+^' (3) 

where 7 is a transfer matrix element. Assuming that only the lowest conduction band 
is occupied by electrons in a semiconducting CNT with u = +1 or —1, we restrict our 
consideration to the n = case at a given temperature T. Then, Eq. (jS)) in the vicinity of 
the K point can be approximated as 

Eg n^e 

with an effective mass m* = l-nfi} j'iL^ and the band gap Eg = 4iT'y/3L. In the K' valley, 
a similar dispersion relation can be obtained when k is substituted by k'. Although it is 
known that the external magnetic field B modifies the CNT electronic states, this effect is 
neglected as the relevant parameter {dt/2aHY (where dt = iLl/n is the CNT diameter and 
an = \/ cfhjeB the magnetic length) is practically very small.i^*ii Hence, we only consider 
the influence of B on electron spin states through the Zeeman energy fiujo] a = ±1/2 is the 
spin projection on the ^ direction. 



Utilizing the expressions given above, we can represent the fluctuating field operator in a 
second-quantized form in terms of the electron creation-annihilation operators aj, ^ and ak,a, 

ki,k2,(7 j=l 

where /i denotes the coordinate for the spin states; by convention, the direction of the 
magnetic field B is chosen as the z axis (quantization axis) and two transversal directions 
as X and y {fi = x,y,z). In addition, a^r = O'hf^o/L and rjj is the location of the j-th 
nuclear spin on the CNT axis. As ki and k2 are any two states in the Brillouin zone, Eq. (0) 
accounts for the effects of both intra- and inter-valley electron scattering on the nuclear 
spins. 

Let us now consider the spin evolution caused by arbitrary random fiuctuations 0^(t). 
The time dependence of the mean spin value s can be described by the quantum kinetic 
equation provided the spin relaxation times Ti and T2 are much longer than the correlation 
time of the thermal bath:— 

j^s{t) = tux sit) -T [sit) - So] , (6) 

where uj = u B /\B \ if the (^-factor anisotropy is ignored and the electron spin polarization 
at thermal equilibrium Sq is given as — |?tanh(/ico'/2A;BT) (A;^ the Boltzmann constant). 
Finally, the matrix T of the relaxation coefficients can be reduced to the Bloch-Redfield 
diagonal form with a leading diagonal composed of matrix elements T^x = ^2"^, ^yy = ^2"^, 
and = Tf ^: 

Tf^ = 27rn(t^)7^^(u;), (7) 
= A^zziO) + niu)-f,xiuj)], (8) 

where nicu) = (1 + e~^f^^'^) /2 and JiM^ico) is the Fourier transformed correlation function 
of the operator 9^, 

= ^J^ (0.(^)0.) e'-^clr. (9) 

Hence, evaluation of the longitudinal Ti and the transversal T2 relaxation times can be 
reduced to finding the relevant 7^^. In Eq. (jHl), ©/^ (t) = exp(ii7rfr/^)9^exp(— zif^T/^), 
(...)= Tr{e~^'^^^'^'^ . . .j/Tre"^'*/*^'^^, where is the Hamiltonian of the thermal bath. In 
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our case, it takes the form 

Hd = ^ ^kflk.a'^k.a + ^ fuJnIz- (10) 
k,CT j 

The first term of Eq. (fTUI) represents the kinetic energy of the electron, which is basically 
the electron Hamiltonian after the Zeeman energy ^cr^k o-'^k.o- is removed; as defined 

k,(T 

earlier, aj^^ and the creation and annihilation operators of an electron with energy 

[Eq- ©] and a is the electron spin quantum number. The second term accounts for the 
magnetic energy due to the nuclear spin splitting in a magnetic field. 

As the electron momentum relaxation time is expected to be shorter than the spin 
relaxation time, the correlation functions can be found from Eq. (jH)) in terms of (5-functions 
reflecting conservation of energy, when the average electron kinetic energy (^k) ~ ksT 
is much larger than the energy broadening T of the order of /^t^^ (i.e., /c^T ^ hr^^). To 
further simply the formulation, the nuclear spin operator Lj. contained in the fluctuating field 
operator 6^; [Eq. (jS}] is conveniently split into two parts Ix = (/++/-) /2 with the raising and 
lowering operators I± = I^^ily] correspondingly, Q± is defined from 6^; = (6+ + 6_)/2 as a 
formal substitution for index fi. Then, by averaging e*'^'^^"^^^''-' over the random distribution 
of nuclear isotopes ^^C, the Fourier transformation 7-1- (uj) of the correlation function 
(e±(r)e^) gives 




The distribution function fk^a = \^(H:,a^k,aj for non-degenerate electrons is 

^-uk,<r/kBT I ^^^^^-Uk,a/kBT ^ whcre Uk^a = Ufc + tvjja, Uk = £k — Eg/2. Since 7++(co') = 

7 (u) = from (/+/+) = {!-!-) = 0, Eq. (fTTj) allows one to find 7xa:(t^) = 

[7_i__(ci;) + 7__|_(co')] /4 as well as n(co')7^.^,(a;) in the form 

n{uj)jxx{^) = ^ [7+-(^) + + 7+-(-^) + 7-+(-^)] • (12) 

Using Eqs. (HH) and ^ and identity {I±I^) = (/J) + (ly) ± {h) = 2 (/J), one can derive 
relaxation parameters in Eqs. ((Tj) and Under the assumption that the nuclear spin 
splitting ujn is negligible compared to c<j, it takes the form 

M^hU^) = i^") [/(I - + /'(I - /)] [Hf^ + e~ e')] , (13) 

k,k',a 
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where e = ek,a, e' = ek',a- Applying inequalities / = fk,a <!,/' = fk',a < 1, Eq- fo^' 
non-degenerate electrons reduces to 

nn{cu)j,,{u) = 2N^ {ll) ^ fk,J{nw + - e^O- (14) 

k,k',u 

Similarly, we find 

vr7..(0) = 2iV^ {II) J2 fk,J{ek - Sk')- (15) 

Note that in the case of / = 1/2, = il that leads to (J^) = (If) = 1/4 (1 is the unity 
matrix) . 

One can see that the contribution of the elastic scattering that does not involve electron 
and nuclear spin fiip-fiop [Eq. ()15p] differs from that of the inelastic process [Eq. ()14|) ] by hw 
in the argument of the 5-function. Subsequently, we focus on the calculation of 7xx(^^) that 
covers the case of Eq. (|15|) in the limit u ^ 0. The sum over the wave vectors in Eqs. (|14p 
and (fT5|) can be calculated by integrating the energy u = Ek — Eg/2 with the density of 
states D{u). In the vicinity of each valley. 



D[u) = ^\ 16 

TTAl V U 



that leads to the electron distribution function in the form 



(7 

It can be shown that the double summation over k and k' can be reduced to the summation 
over a single valley by multiplying the result by the valley degeneracy = 2. Straightforward 
calculation of the integrals with the density of states D{u) under the condition hui <^ ksT 
results in 

/ N / N a^^lvVm* N keT 

In the limit uj —>■ 0, Eq. (fTHj) reveals a logarithmic singularity. This situation is not only 
typical for one- dimensional systems but also known in the galvanomagnetic effect in bulk 
crystals. A standard recipe for removing such a divergency consists of taking into account 
broadening of the energy levels F due to the electron scattering processes discussed earlier. 
Therefore, as soon as hw becomes smaller than this broadening factor F, the magnetic field 
dependence becomes saturated at \ia{kBT/T) instead of In^kBT/fkj) in Eq. (fTHjl . Moreover, 
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one must also take into account the finite length Aq of an actual CNT. In such a case, the 
electron energy cannot be less than Ae ~ h'^/m*AQ, which substitutes F if Ae > T. In 
the following, we assume that the effect of finite Ae is included in the parameter F. Note 
that a similar restriction on the bottom limit of the electron energy u would be applied to 
D{u) in Eq. Therefore the condition for the validity of Eq. (|18|) will be satisfied by 

max{F, fkv} <^ ksT. 

In a manner similar to that discuss above, we can find 72^(0), which looks like Eq. (jlSp 
with hhj — > F. The final expressions for relaxation times [Eqs. ((Tj) and (jH))] take the form 



V = ^(l.^f + (20) 

2 \ F maxjAL;, F| / 

where the essential part, which determines the order of magnitude of the spin relaxation, 
can be expressed in terms of the fundamental CNT parameters 



Equations (fT^ . (pUj) and (PT|) exhibit an unusual temperature dependence for the spin 
relaxation rate; in contrast to the three-dimensional case, decreasing T enhances spin re- 
laxation. Apparently, this effect stems from the property of the one-dimensional density of 
states, which increases as u decreases to F. Equation (PT|) also shows that the geometric 
properties of different CNTs manifest itself only via the length of the chirality vector as a 
factor {L = hi a/^^ + + nin2). Hence, the relaxation rates for a variety of semi- 

conducting CNTs can be readily compared by using this scaling rule. As for the magnetic 
field dependence, it appears in Eqs. (fT^ and (plj) as the parameter f\w, which interplays 
with F. When hw > V, the calculation predicts gradual reduction of the spin relaxation rate 
as B increases. On the other hand, no magnetic field influence can be expected once hjj 
drops below F. 

As an example, we consider a zigzag CNT with {ni,n2) = (8,0) and assume F = 1 fieV. 
Other parameters are known to be: x = 0.011, = -\/3&^/4, b = 0.249 nm, 7 = 70^0/^) 
7o = 3.013 eV.2 The HFI constant for ^^C was estimated in Ref. Q ahf/2'Kh = 22.5 MHz. 
Figure 2 presents the calculated relaxation rates T^^ and Tg"^ as a function of temperature at 
various magnetic field strengths. Clearly, spin relaxation becomes slower with the increasing 
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temperature as discussed above. Ti is always longer than T2 with the exception of the zero- 
field case, where the longitudinal and transversal relaxations are indistinguishable. Both 
relaxation rates also show gradual decrease as B becomes larger. With the relaxation time 
of about 1 s, these characteristics are readily observable by experiments. 

In order to consider the effect of radial confinement, we calculate the spin-fiip rate Wd in 
bulk diamond. In general, Wd = nvcTsf, where n = 2x/fld is the nuclear spin concentration 
{Qd is the unit cell volume of diamond), v = ^/^k^TJrnri the mean electron velocity at 
temperature T, and m = ^/m^jn^ the density of states effective mass, m,, and m± are the 
longitudinal and transversal masses in each of diamond X valleys with six fold degeneracy 
(i.e., ly = 6). The spin-flip cross section for electron scattering with a localized spin moment 
calculated in the flrst Born approximation is known to be asf = -^^1(1 + l)/^a|jf2^m^ [see 
Ref . Iisl . Taking into account that Qd = 5.67 x 10~^^ cm^, m,, = IAitlq and m± = 0.36mo (mo 
is the free electron mass), one can estimate the spin-flip relaxation time W^^ = 4.7 x 10^ s. 
This value exceeds the CNT relaxation time at T = 4 K and S = by at least of four orders 
of magnitude, demonstrating the signiflcance of the radial conflnement effect in a CNT. 

In conclusion, we consider electron spin relaxation in a single-walled semiconducting CNT 
through the HFI with nuclear spins of ^^C isotopes. The analysis reveals the peculiarities 
in spin relaxation inherent to one dimensional systems at low temperatures and/or weak 
magnetic flelds. As a result, it becomes dependent on the non-magnetic electron scattering. 
Numerical estimations illustrate the relative importance of this relaxation mechanism in 
a CNT compared to the similar processes in bulk diamond and other carbon-based struc- 
tures; strong enhancement due to the radial conflnement of electrons helps making the HFI 
dominant over the spin-orbital interactions, particularly at weak magnetic flelds and low 
temperatures. However, the anticipated spin relaxation time of the order of 1 s in CNTs is 
still much longer than those found in conventional semiconductor structures. 

This work was supported in part by the SRC/MARCO Center on FENA and US Army 
Research Office. 
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FIG. 1: Upper: Lattice structure of the graphene sheet. The carbon atoms are located at the 
vertices of hexahedrons, b i and b 2 are primitive translation vectors. L is the chiral vector 
and denotes the direction perpendicular to L . The lengths of the vectors are b = \/3ac-c] 
L = h\lri\ + + n\n2\ ac-C is the distance between the nearest carbon atoms. The figure depicts 
the particular case of ni = 4, 71-2 = 2. Lower: CNT as a tortile graphene sheet. ^ and 77 denote the 
coordinates for the electronic states. The direction of the magnetic field ^ is also shown. 
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FIG. 2: Calculated spin relaxation rates (a) T[' and (b) T2 in a (8,0) zigzag CNT as a function of 
temperature for different values of magnetic field B. The strength of the magnetic field is indicated 
in units of Tesla. 
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